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AN ANALOG STUDY OF THE PASSIVE THERMAL BEHAVIOR 
OF AN ORBITING SPACE STATION 

By Lenwood G. Clark and John W. Wilson 
Langley Research Center 

SUMMARY 

A relatively simple procedure for determining the passive thermal-control capa- 
bility of an example space station is presented. Approximate Fourier heat equations 
governing the heat flow through the assumed multilayer wall structure of the spacecraft 
were solved by means of a general-purpose analog computer. The thermal behavior of 
the spacecraft was determined for a number of surface -coating distributions. Of these, 
the surface -coating distribution which provided the least circumferential temperature 
gradient was selected and the passive thermal-control capability of the spacecraft was 
determined for a range of orbit conditions. Comparison of steady-state and transient 
solutions showed that the steady-state solution underestimated the passive heat-transfer 
capability of the spacecraft. 

The analog computer, with its usual ease of changing parameters and rapidity of 
solution, proved to be relatively economical. An analytical solution of linearized equa- 
tions was in good agreement with the corresponding analog solution. 

INTRODUCTION 

For the past decade, increasing amounts of the literature have indicated the impor- 
tance of temperature in the design of a spacecraft. The thermal design of any spacecraft 
is subject to the stringent temperature limitations imposed by each of its onboard compo- 
nents, whether they be man or machine. The degree to which the thermal design is inte- 
grated into the structural design depends not only upon the required component tempera- 
ture regimes, but also upon the method of thermal control used. 

From the reliability standpoint, passive thermal control achieved at the time of 
construction by the application of proper optical coatings, careful selection of the space- 
craft geometry, and judicious selection of construction materials (ref. 1) would be desir- 
able. However, in many instances it must be supplemented by active thermal-control 
devices to achieve the required temperature levels. 



The purpose of this report is to present a relatively simple method by which esti- 
mates can be made of the passive thermal- control capability for providing an optimum 
living and working thermal environment for the crew members of a manned spacecraft. 

As an example, the report considers the thermal control during the long-term orbital 
phase (1 to 5 years) of a large (150-foot-diameter (45.7-m)) hexagonal space station 
concept shown in figure 1. Knowledge of the thermal mission requirements of the vehicle 
will permit estimates of the degree of active thermal control which will be required to 
maintain desired temperatures. On this basis, the vehicle can either be redesigned, if it 
appears necessary, or the present design can be analyzed more rigorously to define its 
thermal-control requirements more precisely. 

SYMBOLS 

The units for the physical quantities used in this paper are given both in U.S. 
Customary Units and in the International System (SI). Factors relating these two systems 
of units are given in reference 2. 

A semimajor axis of orbit ellipse, nautical miles 

a earth albedo, dimensionless 

C thermal capacitance per unit area, Axpc, Btu/ft2-°F (j/m 2-°c) 

C s solar radiation constant, Btu/ft2-hr (w/m2) 

c specific heat, Btu/lbm-°F (j/kg-°C) 

D vector from spacecraft to center of earth 

e orbit eccentricity, dimensionless 

F geometric shape factor for a spinning sun-oriented planar element, 

dimensionless 

H distance between spacecraft and center of earth, nautical miles 

h orbit altitude, nautical miles 

h c forced convection heat -transfer coefficient, Btu/4t2-hr-°F fj/m^-s-°c) 
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I orbit inclination measured from equatorial plane, deg 

i ecliptic inclination measured from equatorial plane, deg 

K earth gravitational field constant, nautical miles^/hr^ 

k thermal conductivity, Btu/ft-hr-°F (j/m-s-°C) 

L(£) length of cylindrical module at circumferential position | , ft (m) 

l wall thickness, ft (m) 

M spacecraft angular-momentum vector 

N(£) vector normal to spacecraft surface element at circumferential position £ 

Q radiant heat flux absorbed by spacecraft surface area, Btu/ft^-hr (w/m^) 

Q tot total radiant heat flux absorbed by spacecraft, Q X Area, Btu/hr (W) 

q radiant heat flux incident on spacecraft surface area, Btu/ft^-hr (w/m^ ) 

q c passive heat transfer rate, Btu/hr (W) 

R thermal resistance, Ax/k, ft2-hr-°F/Btu (m^-s- °c/j) 

Rp; mean radius of earth, nautical miles 

r radius, ft (m) 

T temperature, °F (°C) 

T temperature, °R (°K) 

T time derivative of temperature, °F/hr (°C/s) 

t time, hr (s) 

X,Y,Z inertial coordinate frame where Z-axis coincides with north pole and 35-axis 
is alined with vernal equinox (see fig. 2) 
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X',Y',Z’ rectangular coordinate frame with Z'-axis normal to orbit plane and X'-axis 
coinciding with earth-spacecraft line (see fig. 2) 

X",Y",Z" rectangular coordinate frame with Z"-axis normal to ecliptic plane and 
X"-axis coinciding with earth -sun line (see fig. 2) 

x wall coordinate, ft (m) 

a absorptivity to solar radiation, dimensionless 

a L angle between line from spacecraft to center of earth and line from space- 

craft to element on earth surface at terminator, deg 

(3 right ascension of ascending orbit node measured in equatorial plane (zero 

at vernal equinox), deg 

y spacecraft position angle measured in orbit plane (zero at ascending orbit 

node), y p + 77, deg 

7p argument of orbit perigee measured in orbit plane (zero at ascending orbit 

node), deg 

A finite difference operator, dimensionless 

e total hemispherical emissivity, dimensionless 

77 spacecraft position angle measured in orbit plane (zero at orbit perigee), deg 

77_ angle between line from center of earth to sun and line from center of earth 

to spacecraft, deg 

| angle between vector normal to spacecraft surface element and spacecraft 

angular momentum vector, deg 

p density, lbm/ft^ (kg/m3) 

c r Stephan -Boltzmann radiation constant, Btu/ft 2 -hr-°R 4 (w/m 2 -°K 4 ) 

a ) angular velocity, radians/sec 
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Q, position angle of sun measured in the ecliptic plane (zero at vernal equinox), 

deg 

right-hand limit as x approaches x 0 
left-hand limit as x approaches Xq 


Subscripts: 


A l 

aluminum 

air 

cabin air 

e 

earth thermal 

ent 

shadow entrance 

exit 

shadow exit 

f 

foam 

i 

inner wall surface 

j 

integer 

N,n 

integers 

P 

perigee 

r 

earth-reflected solar 

s 

solar 

space 

space environment 




ANALYSIS 


Development of Heat -Flow Equations 

The space station consists of six cylindrical modules arranged in the shape of a 
hexagon and connected by three spokes to a central hub which provides docking facilities. 
(See fig. 1 and ref. 3.) The normal mode of operation calls for the station to be spinning 
about a vector normal to the plane of the hexagon and alined toward the sun. Because of 
symmetry in the hexagonal station and because negligible direct and reflected radiation 
from one portion of the station to another is assumed (since subtended angles are gener- 
ally small), one of the cylindrical modules is representative of the complete station. 

A schematic drawing of the assumed basic wall structure (excluding structural 
members such as ring stiffeners and longerons) is shown in figure 3. As shown, the 
wall is essentially a five-layer composite. On the outside is a thick layer of polyure- 
thane foam; on the inside is a thick layer of aluminum honeycomb with a polyurethane- 
foam-filled core. These layers are separated and covered inside and out by thin 
aluminum sheets. 


During an equilibrium orbit, the net heat flow between the internal environment of 
a module and the enclosing wall surface equals the net heat flow through its wall struc- 
ture. If the former heat transfer is assumed to occur solely by forced convection (this 
assumption is discussed later in connection with boundary conditions), the latter (passive) 
heat -transfer rate through the wall is given by 


% 



Ti(*,t) 


d£ dt 


( 1 ) 


where from figure 4 the length L(£) of the cylindrical module at circumferential posi- 
tion £ is given, 


in feet, by L(£) = 69.2269 + 5.7731 sin £ 


in meters, by L(£) = 21.1004 + 1.7596 sin i; 


( 2 ) 


and the inner wall surface temperature T^(£,t) is found by solving Fourier's heat equa- 
tion subject to appropriate boundary conditions. 

If negligible heat flow is assumed in the circumferential and longitudinal directions 
and the slight wall curvature is neglected, the Fourier heat equation reduces to 

8 2 T _ pc 8T /„ 

q x 2 k 8t ' 
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In terms of the physical wall illustrated in figure 5(a), equation (3) requires that there be 
two paths of heat conduction through the honeycomb with no cross conduction between the 
two. 

If a digital computer is employed in the study, both independent variables must be 
differenced. However, since an analog computer was to be used to solve the equations, 
there existed an option for allowing one of the independent variables to be continuous. In 
this study, the time variable was made continuous and finite differences were used for the 
linear dimension. This choice not only assures the numerical stability of the solution, 
but allows convenient handling of the boundary conditions (ref. 4) and the thermal proper- 
ties of the wall which are discontinuous functions of position. 

If a second-order correct central difference approximation to the spatial derivative 
is employed, equation (3) becomes 


gc 9T 
k 0t 


= ( T j-1 - 2 Tj + Tj +1 ) / ? \ 

8x2 lx. (A Xj )2 ^ l) 


(4) 


If the remainders 0 (ax?^ are neglected, the resulting differential equations governing 
the heat flow per unit area through the wall in a radial direction for the nine selected 
node positions of figure 5(b) can be written as 


V^fr-i- 2T i + T i+i) < 5 > 

The values of thermal resistance and thermal capacitance C n are shown in 
table II; these values were computed from the thermal properties (assumed constant for 
this simplified analysis) and physical dimensions of the wall presented in table I. The 
increment sizes must be chosen to bound the truncated remainder to a small value. It is 
interesting to note that the solutions on an analog computer with t continuous appear 
"satisfactory even for rather large values of Ax" (ref. 5) as opposed to a digital- 
computer solution where the increment sizes ^At and Axjj must be chosen not only to 
bound the two truncation errors, but in some cases to insure numerical stability. 

Each layer of the wall structure is thermally coupled to its environment. This 
thermal coupling is expressed as a set of boundary conditions describing the mode of 
heat transfer at the interior-wall surface, the interfaces of the composite layers, and the 
exterior-wall surface. Consideration of these boundary conditions yields additional 
equations which can be approximated by first-order correct forward and backward 
differences. 
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Interior wall. - If all internal heat transfer (i.e., metabolic and equipment heat) to 
the interior wall surface is assumed to occur by convection, the boundary condition to be 
satisfied is 

kM £U_> - T «) + 0(ax) = - T is) < 6 > 

This assumption was made because of the many uncertainties regarding the time- 
dependent internal distribution of onboard heat -producing components; which, even if 
known, would be difficult to describe. However, consideration of human comfort will 
undoubtedly influence the dependence on forced convection as a mode of heat transfer in 
manned spacecraft. 

Interfaces .- If zero contact resistance is assumed at all interfaces between mate- 
rials, the boundary condition to be satisfied at the jth interface is then 

k . »T| =k . . 9T 

j 8 x lx=Xj(-) 9x1 X=Xj(+) 

where the subscripts j and j+1 refer to the two dissimilar materials forming the 
interface. This assumption is reasonable if it is assumed that a good bond is present at 
these interfaces. Realistically, the joint resistances at these interfaces should be deter- 
mined experimentally. Unfortunately, the mode of heat transfer through these interfaces 
is not well understood and prior experimental data regarding this problem have proved 
unreliable because of data scatter and the general inability to reproduce experimental 
results (ref. 6). The resulting boundary condition obtained from finite difference approxi- 
mation is 



Exterior wall .- When all heat transfer at the exterior wall surface is assumed to 
occur by radiation, the boundary condition to be satisfied is 

Com S ' Tl ) = Q({,t) ■ w ( T * ■ d <8) 

where Tf„ 0 „„ is small and can be neglected, and Q(£,t), the absorption rate of 

Sp3.Cc 

external radiant energy from the sun, the earth, and sunlight reflected by the earth, is 
discussed in appendix A. 
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To aid in the solution of equations (1) and (5), the cylindrical module is approxi- 
mated by a 12 -sided polyhedron as shown in figure 4; equation (1) can then be approxi- 
mated by 

12 

<lc = £ <1 c(Sn) ( 9 > 

N=1 

where 

1c(*N) = ^*{f)L(«N) £ fair - T 15 (« N ,t)|dt (10) 

For the spinning case, the incident external thermal energy is symmetric about the 
angular -momentum vector and only six sides need be considered (i.e., = 15°, 45°, 

75°, 105°, 135°, and 165°). 


Analog-Computer Equations 

Equation (5) for the first aluminum sheet (j = 2, n = 1) places serious burdens on 
the computer used, analog or digital. The coefficient (l/RjCi « 10 6 /hr) is large, and 
the temperature differences are small, that is, [JTj - 2T2 + TgJ * 10" 5 °F or °c]. 

Hence, much precision is required to evaluate the derivative without loss of significance. 
Perturbation about some mean temperature, although helpful, will not eliminate this diffi- 
culty. Time scaling the analog computer is of no practical use because of the great dis- 
similarities in thermal properties and linear dimension of the various wall layers as 
shown in tables I and n. 

A general technique to eliminate this loss of significance is described in appendix B. 
The technique involves judicious substitution of the boundary conditions into equation (5), 
after which the following simplifying assumptions can be made without introducing serious 
errors into the computation: 


T 1 = t 2 = t 3 
t 6 = t 7 


Tl4 = T 15 


If the remainders are neglected, the computer equations (eqs. (5) to (8) and eq. (10)) are 

Tl = K t [K 2 Q(? N ,t) - ae t| + K 3 (t 4 - T^j (11) 

T 4 =K 4 (t 1 -2T 4 + T 5 ) (12) 
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( 13 ) 


T 5 = K 4 (T 4 - 2T 5 + T 6 ) 

T 6 = K 5 [k 6 T 5 - (2 - K ? - K 8 )T 6 + K 9 T 9 + K 10 T^ (14) 

T 8 =K 8 T 6 + K 9 T 9+ K 10 T 11 (15) 

T 9 =Kii(t 8 -2T9 + T 10 ) ( 16 ) 

ti 0 = K 11 (T 9 -2T 10 + T 13 ) (17) 

Tll=K 12 (T 8 -2T 11+ T 12 ) (18) 

fi2=Ki 2 (T 11 -2T 12+ T 13 ) (19) 

Tis = K 8 T 15 + K 9 T 10 + K 10 T 12 (20) 

*15 = K 5 [k 9 T 10 + K 10 T 12 - (2 - Kg - K 13 )T 15 + K 14 T ai ^j (21) 

and t 

9c(In) = k 15(«n) J 0 |_ T air - T 15(«N.‘J d ‘ < 22 > 


where the coefficients K n are defined in terms of the equivalent network parameters by 


Kl 


_ 1 
R l c l 


K 2 = R 1 


K 3 = 


R 1 

R 1 + R 2 


k 4 = 


1 

r 2 c 2 


k 5 


1 

r 3 c 3 


r 6 


R 2 

r 2 + r 3 


k 7 = 


r 3 

R2 + r 3 


k 8 = 


r 4 r 5 

R^Rg + R3R5 + R3R4 


k 9 = 


r 3 r 5 

R4R5 + R3R5 + R3R4 


Kiq = 


r 3 r 4 

R^Rg + RgRg + R3R4 


Kll = 


1 

R 4 C 4 


Kl2 


1 

r 5 c 2 


K13 


r 3 

r 3 + r 6 


Kl4 


R6 

r 3 +r 6 
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and 


Equations (11) to (22) were programed on a general-purpose analog computer. The 
analog circuit for this system of equations is shown schematically in figure 6 where it 
should be noted that the coefficients of T-q and in equations (15) and (20) were 

too small for consideration. The absorption rate of external thermal energy Q(^,t) 
was supplied by a function generator which was programed from a digital-computer solu- 
tion of equations (Al) to (A10). 

To verify the analog circuitry, the governing equations were linearized by replacing 
the term Q^N>i) " oet| by Q cos cot ^where Q = 0.2778 Btu/ft^ -sec ^3150 j/m^-sec) 
and co = 0.0714 rad/sec). The linear set of equations was then solved by using the 
analog circuit shown in figure 6 (by disconnecting the multipliers and programing 
Q cos cot on the function generator). These equations were also solved analytically by 
using the impulse response technique (one-dimensional Green's function). See 
reference 7. 

Comparison of the two solutions is shown in figure 7. The small difference in the 
solutions is attributed to slight errors in that part of the analog circuit which excluded 
the multipliers. 


RESULTS AND DISCUSSION 

The system of equations (11) to (22) derived for a planar element of the composite 
wall was applied to each of the planar sides positioned by ijjj for several surface optical 
property combinations. Figure 8 shows the resultant outer-wall surface-temperature 
histories for a circular 200-nautical-mile minimum sunlight orbit (orbit wherein the posi- 
tion of the sun is colinear with the intersection of the orbital and ecliptic planes) and a 
constant air temperature of 70° F (21° C). As illustrated by the figure, a wide variation 
in temperature is possible for each side depending on the optical coating selected. 

If, at a given instant, a large circumferential temperature gradient exists, the 
resultant thermal stress may cause delamination of the bonded joint between wall layers. 

If such a delamination occurred, it would not only alter the thermal characteristics of the 
wall, but might affect the structural integrity of the vehicle. Thus, a coating distribution 
was selected to reduce the circumferential temperature gradient. This distribution is 
shown in table in. Figure 9(a) shows the outer-wall surface -temperature distribution for 
the selected set of coatings. The corresponding internal-wall surface -temperature dis- 
tribution is shown in figure 9(b). 
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The noon orbital position (point of the closest approach of the vehicle to the sun) in 
figure 9 was selected as representative of an extreme circumferential temperature gra- 
dient in order to examine the temperature distribution through the wall structure of the 
vehicle. Figure 10 shows this result. Of interest is the small temperature differences 
which exist across the face of the thin aluminum sheets. In an a priori fashion this fact 
justifies the previous assumption of small temperature differences used in developing 
the equations. 

For the selected coating distribution, the external and internal surface temperature 
histories were then calculated for a circular 300-nautical-mile orbit for both minimum 
and maximum sunlight (orbit wherein the position of the sun is normal to the intersection 
of the orbital and ecliptic planes) conditions to bracket the design range for this space 
station completely. These results are shown in figures 11 and 12 for a constant air tem- 
perature of 70° F (210 C ). 

The passive heat -transfer capability of the cylindrical module with the selected 
optical coating distribution is shown as a function of the heat absorbed by the module in 
figure 13. Of interest, is the steady-state solution (based on an average external-wall 
surface temperature) which is shown for comparison. The figure shows that a steady- 
state solution underestimates the passive heat -transfer capability of the vehicle. 

Figure 14 shows the heat absorbed by the cylindrical module having the selected 
qptical coating distribution as a function of altitude for both minimum and maximum sun- 
light orbits. This figure indicates that for a given altitude the vehicle absorbs less heat 
for a maximum sunlight orbit than for a minimum sunlight orbit for the vehicle orienta- 
tion considered. The reason for this absorption difference is that the vehicle is oriented 
during a maximum sunlight orbit so that it receives less earth-thermal and earth- 
reflected solar radiation. Together, figures 13 and 14 define the passive thermal-control 
capability of the cylindrical module (with the selected optical coating distribution) for the 
altitude range of interest. 

It should be noted that equations (11) to (22) are in a form to be integrated by using 
standard subroutines on a digital computer. However, past experience has shown that 
such practice is expensive because of the long computation time which is due to the use 
of extremely accurate integration subroutines with too stringent error criteria. If a 
digital solution is considered, then a simpler approximation to the time derivative is 
suggested by reference 4. 


CONCLUDING REMARKS 

The passive thermal behavior of an example space station was determined for a 
range of orbit conditions from analog-computer solutions of Fourier's heat equation as 
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applied to the composite multilayer wall structure of the spacecraft. Comparison of 
steady-state and transient solutions showed that the steady-state solution underestimated 
the passive heat-transfer capability of the spacecraft. 

Explicit use was made of Fourier's heat equation which reduced the chance of an 
erroneous solution that may result from intuitive approximations. Approximations such 
as lumping the aluminum with the foam, neglecting the thermal conductivity of the alumi- 
num, and even neglecting completely the aluminum sheets, are not necessary and give no 
indication of the amount of error introduced. In this report, the errors other than trunca- 
tion errors, which are predictable, are known to be compatible with the resolution of the 
general-purpose analog computer. 

The governing heat -transfer equation can be programed directly if the composite 
wall layers are similar in both linear dimension and thermal properties. When great 
dissimilarities exist, loss of significance will occur. This loss of significance can be 
relieved by the technique outlined in an appendix. 

The analog computer, with its usual ease of changing parameters and rapidity of 
solution, proved to be relatively economical. An analytical solution of linearized equa- 
tions was in good agreement with the corresponding analog solution. 

Results of this study are of value in the thermal design of any orbiting spacecraft 
since they illustrate a relatively simple procedure for determining the passive heat- 
transfer capability of a vehicle. Knowledge of the thermal mission requirements of the 
vehicle will then permit estimates of the supplemental active thermal control necessary 
to maintain desired temperatures. On this basis, the vehicle can either be redesigned if 
it appears necessary, or the present design can be analyzed more rigorously to define its 
thermal control requirements more precisely. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., June 17, 1966. 
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APPENDIX A 


DETERMINATION OF EXTERNAL THERMAL ENERGY 


The external thermal energy which is absorbed by a spacecraft surface element is 
given as a function of element position £ and orbit position 77 by 


QU,V) = [%(!,??) + q r (|,77)]o! + (Al) 

where 

q = C s cos | (|l|<90°); (n exit <T]< 

<l s (Z>v) = 0 (Nl =90°); (?7 ent ^77 ^ 



and 


q r (£, 7 ?) S C g aF( 4 , 77 )cos t] s 
q r (^?7) = 0 



4 


(A4) 


Time -Position Relation 


The absorbed external energy Q(£,? 7 ) as defined by equation (Al) can be expressed 
as a function of time in orbit (measured from orbit perigee) to reach orbit position 77 
(i.e., Q(|,t)). This relation is given from reference 8 for an elliptical orbit by 


t = 



1 /e + cos 77 | ^ _ eyl - e2 sin 77 

1 1 + e cos 77 / 2 1 + e cos 77 


(A5) 


Geometric Shape Factor 

From reference 9, the geometric shape factor F(£,7j) for a spinning sun-oriented 
planar element in equations (A3) and (A4) is given by 
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APPENDIX A 


F(4,t?) = B 0 + B^C) + B 2 (c 2 + ^ + Bg ^C 4 + 3C2 d 2 + |d 4 


b 4 c 


+ i£c 4 D 2 + ^C 2 D 4 + 

2 o 16 J 


(A 6 ) 


where 


C = -cos | cos rj s 
D = sin | sin r 7 s 


and where 


> 

B ° = |^(l05 " 7 cos “L + 3 cos 3 Q! L - | cos5a L + 1 cos 7 ^ 


Bi = | sin 2 ^ 


Bg = t^( cos a ij ~ 2 cos 3 q!l + 4 cos^ol - 3 cos^aj^ 
b 3 = t^(-cos “l + cos 3 a L - ^ cos 5 Q! L + 18 cos 7 ®^ 
B 4 = gg^ (p cos q-l - 35 cos 3 a L + 63 cos 5 q! L - 33 cos 7 ^ 


• -i /R] 
= sin 1 — i 


0 !l = sin 


V« 


! ) 


H = 


_ ( h p + r e )( 4 + e ) 


1 + e cos it 


h = H - R E 


(A7) 


Attitude -Position Relation 

The attitude angle 180° - r\ s shown in figure 2 defined in terms of orbit param- 
eters for a sun-oriented vehicle can be found since (see ref. 9 for derivation) 


15 



APPENDIX A 


cos 7 j s = cos S2(cos y cos /3 - sin y cos I sin 0) 

+ (cos i sin S2)(cos y sin /3 + sin y cos I cos j3) 

+ sin i sin sin y sin I (A 8 ) 

where it should be noted that y = + 77 . 

Shadow Entrance and Exit 
Shadow entrance and exit occur when (from ref. 10): 


cos r] s = -cos Q!l 

where 



(A9) 


(A10) 


The limits on the position of the spacecraft for receiving both solar radiation (eqs. (A2)) 
and earth-reflected solar radiation (eqs. (A3)) are then determined from equations (A 8 ), 
(A9), and (A10). 
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APPENDIX B 


TECHNIQUE FOR INCREASING RESOLUTION IN COMPUTATIONS 


Differential equation (5) for n = 1, 3 and j = 2, 5, 14 govern the heat flow in the 
thin, highly conductive aluminum layers of the wall structure. Significance can be 
increased by judicious substitution of the algebraic boundary conditions (eqs. ( 6 ), (7), 
and ( 8 )). To illustrate this, consider the following forms of equations (5), (7), and ( 8 ): 

t 2 = kA * ro( T l - 2 T 2 + t 3 ) 

P Al c Al( Ax l) 

(Bl) 

^( T 3- T a)=^j( T 4- T 3) 

(B2) 

- t 2 ) = Q(£,« - 

(B3) 


where Tj and T 3 are the temperatures on the boundary and T 2 and T 4 are the 
temperatures inside the materials and are described by time derivatives. The conditions 
of dissimilarity in thermal properties and physical dimensions are 

kA l » kf A Xl « Ax 2 

— » 1 ^Al » 1 

pai c ai( Ax i) Ax ! 


Let bi = bo = and ai = — - — so that bj » b 2 , aj » 1, and 

AX 1 AX 2 (Ax^PAfC^ 

bj » 1. 

It follows that 


i « 1 ^2 « 1 (B4) 

b l b l 

Solving equations (B2) and (B3) for T 2 and Tj, respectively, and substituting into equa- 
tion (Bl) gives 
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APPENDIX B 


T 2 = j^Q(S,t) - oeT\ + b 2 (T 4 - Tg)j (B5) 

where the coefficients for this equation are ^ and From the inequalities 

b l b l 

(eq. (B4)), it follows that ^ « a^ and « aj. The coefficients in the new form of 

equation (3), equation (B5), have thus been greatly reduced and the differences in the 
bracket are larger. 
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TABLE L- MATERIAL PROPERTIES 



p 

k 


c 


Material 

lbm 

kg 

Btu 

J 

Btu 

J 


ft3 

m3 

hr-ft-°F 

m-s-°C 

lbm-°F 

kg-°C 

Aluminum 

173 

2770 

110 

190 

0.211 

833 

Polyurethane foam 

1.50 

24.0 

.009 

.016 

.250 

1050 

Aluminum honeycomb 

4.50 

72.1 



— 

— 


TABLE n.- HEAT -TRANSFER CHARACTERISTICS 
AND PHYSICAL DIMENSIONS 



Ax n 

Rn 

c n 

n 

ft 

m 

ft2 -hr-°F 

m 2-s-°C 

Btu 

J 


Btu 

J 

ft 2 -°F 

m2-°C 

i 

8.33 X 10" 4 

2.54 X lO" 4 

7.57 x lO" 6 

1.34 X 10-6 

3.04 x 10-2 

5.86 x 10 3 

2 

2.64 x lO- 2 

8.05 X lO" 3 

2.93 

5.03 X 10" 1 

9.90 x 10-3 

2.02 X 103 

3 

1.58 x lO” 3 

4.82 x 10-4 

1.44 X 10-5 

2.54 x 10-6 

5.77 X 10-2 

1.11 X 104 

4 



9.23 x 10-3 

1.63 X lO" 3 

9.58 x 10' 1 

1.86 x 10 5 

5 



3.01 

5.16 x lO" 1 



6 



1.00 

1.76 x 10" 1 




TABLE HI.- SELECTED THERMAL COATING DISTRIBUTION 


£n 

15° 

450 

75° 

105° 

135° 

165° 

O' 

0.3 

0.3 

0.5 

0.5 

0.3 

0.3 

€ 

.9 

.9 

.5 

.1 

.1 

.1 
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Figure L- Artist's concept of a 150-foot-dia meter (45.7-m) hexagonal space station. L-62-8400 
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(a) Physical system shewing temperature node locations. 
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(b) Basic electrothermal circuit. 

Figure 5.- Schematic of composite thermal wall, 
to 

CJI 




Figure 6.- Schematic of analog-computer circuit for determining the thermal behavior of a wall element of the cylindrical module. 
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Figure 8.- Continued. 
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(a) Outer wall. 
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(b) Inner wall. 


Figure 11.- Temperature history and temperature distribution for minimum sun as a function of orbit position for selected thermal coatings. 

(See table III.) 
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Figure 13.- Passive heat transfer as a function of heat absorbed for cylindrical module with selected coating distribution. (See table III.) 
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